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ABSTRACT 

Context. Fairall 51 is a polar-scattered Seyfert 1 galaxy, a type of active galaxies believed to represent a bridge between unobscured 
type-1 and obscured type-2 objects. Fairall 51 has shown complex and variable X-ray absorption but only little is known about its 
origin. 

Aims. In our research, we observed Fairall 51 with the X-ray satellite Suzaku in order to constrain a characteristic time-scale of its 
variability. 

Methods. We performed timing and spectral analysis of four observations separated by 1.5, 2 and 5.5 day intervals. 

Results. We found that the 0.5-50 keV broadband X-ray spectra are dominated by a primary power-law emission (with the photon 
index « 2). This emission is affected by at least three absorbers with different ionisations (log£ » 1—43- The spectrum is further shaped 
by a reprocessed emission, possibly coming from two regions - the accretion disc and a more distant scattering region. The accretion 
disc emission is smeared by the relativistic effects, from which we measured the spin of the black hole as u * 0.8 ± 0.2. We found 
that most of the spectral variability can be attributed to the least ionised absorber whose column density changed by a factor of two 
between the first (highest-flux) and the last (lowest-flux) observation. 

Conclusions. A week-long scale of the variability indicates that the absorber is located at the distance 0.05 pc from the centre, i.e., 
in the Broad Fine Region. 
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1. Introduction 

In the unification paradigm of active galactic nuclei (AGN) by 
Antonucci ( 1993 !). type 1 and 2 AGN are distinguished by the 
inclination angle, more specifically whether or not the torus in¬ 
tercepts the line-of-sight. However, this standard picture assum¬ 
ing a h omogeneous dusty torus has been questioned by severa l 
studies (iNenkova et al J[2QQ2klElitzurll2Q 1 2llMerloni et al.ll2013h . 
The clumpiness of the absorber is supposed to be another im¬ 
portant parameter for the proper classification. The type of AGN 
is not entirely determined by the inclination angle but also by 
a probability of the absorption clouds intercepting the line of 
sight. This is supported by a large and fast variability of the_ab- 
sorber observed in s everal sources (see, e.g jRisaliti et al.l[2005b 
Bian chi et al.ll2012h . The statistica l analysis based on the vast 
RXTE archive was recently done by Marko witz et al.1 (12014l) giv¬ 
ing the probability of an absorption event regardless the con¬ 
stant absorption due to non-clumpy material to be 0.006'[[' 160 


.003 


for type 1, and 0.110^ ^j for type 2 AGN, respectively. 

The absorption variability occurs at different time scales 
for different AGN. iRisaliti et al.l ( 2002 1 investigated a sample 
of Seyfert 2 galaxies and revealed that most of them are 
variable on time scales of months to years. Similar variability 
was recently reported based on Suzaku observations, e.g., by 
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iBraito et al.l (l2013h for NGC 4507, and iRivers et al.l (12014 
for NGC 2110. The variability time scale of months to years 
is typical for a clumpy torus at the distance ~ lpc. However, 
the absorption variability was later found on the time-scale 


of days in several sourc e s: NG C 4388 (lElvis et al.l 20041), 


NGC 1365 (IRisaliti et al.l 120051). NGC 41 51 dPuccetti et al.l 


l2007h. NGC 75 82 (iBianchi et al.l l2009bl). PG 1535+547 
| Balloetal.ll2008h , Mrk 766 (IRisaliti etalJ 1201 lh. NG C 5506 

— ET" 


Markowitz etal.1 I20l4 . NGC 3227 (La mer et al.l 120031 


Beuchert et al. 2015), or Swift J2127.4+5654 (Sanfrutos et al 


20131) . The origin of such varia bility must be significantly closer 
to the central region. Recently. iTorricelli-Ciamponi et al.l ( 2014 ) 
showed that X-ray eclipses are common when the expected 
occultation time is compatible with the exposure and that these 
events may be responsible for most of the spectral variability 
at energies higher than 2 keV. They associated the origin of the 
eclipsing clouds to the Broad Line Region of AGN. 


Knowledge of the location of the absorber in the so called 
polar-scattered Seyfert 1 AGN is important to understand the na¬ 
ture of this kind of objects. Generally, Seyfert 1 galaxies do not 
exhibit polarised emission. It is believed that the total polarisa¬ 
tion is washed out by a mutual interaction of polarised emission 
from two s cattering regions - the equatorial plane an d the ionisa¬ 
tion cone (ISmith et al.1 2004t iBatcheldor et al.ll201 lh . The equa¬ 
torial emission is obscured by a dusty torus in Seyfert 2 galaxies 
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and thus, the polarised e missio n is detected only from the ioni¬ 
sation cone. ISmith et ali (120041) argued that the detected polari¬ 
sation is inclination-dependent and that polar-scattered Seyfert 1 
galaxies represent a bridge between type 1 and type 2 galaxies. 
In this scenario, the enhanced variability of the absorber (with 
respect to the equatorial-scattered ones) is therefore expected be¬ 
cause we may be seeing the AGN along a line-of-sight grazing 
the torus rim. 

X-ray properties of polar-scattered Seyfert 1 galaxies were 
studied bv lJimenez-Bailon et al.1 ( 2008 1 who realised that the lu¬ 
minosity and the index of a power-law emission are consistent 
with Seyfert 1 galaxies. A low inclination is also suggested from 
detection of a relativistic ally smeared iron line profile coming 
from the innermost accretion disc in most of the sources. The 
variability and complex absorption are common features of these 
sources. 

More recently, investigation of the variability properties was 
done for several sources. A clear example of the absorption vari¬ 
ability was show n in the case of Mrk 704 dMatt et alJ 1201 lb 
lLaha et alil201 11 1 where the spectra from two different obser¬ 
vations separated by three years were identical above 7 keV 
but very different at lower energies. The observed X-ray spec¬ 
tral variability of Mrk 231 be low 10 keV can be interpreted as 
changes in a patchy absorbe r dPiconcelli et al.ll201 3b [Teng et al. 
120141) . Beuche rt et alJ (1201 5l) found the presence of a non-dusty 
absorber in NGC 3227 varying on a time scale of about one 
week. They found that the absorber is intermediately ionised 
and partially covering the source. The spectral variability of an¬ 
other polar-scattered Seyfert 1 galaxy, ESO 323-G77, has been 
recently studied bv iMiniutti et all ( 2014 1. They found two vari¬ 
able absorbers with a different time-scale of their variability. The 
absorbers were associated to Broad-Line-Region clouds (month¬ 
long time scale) and a clumpy torus (year-long time scale) ac¬ 
cording to their variability. 

Fairall 51 (also known as ESO 140-43) is a nearby galaxy 
with the cosmological redshift z = 0.0141, as measured in the 
2MASS survey dHuchra et al. 2012 ). Fairall 51 is classified as 
Seyfert 1 since its optical spectrum contains broad lines with 
the full width at half maximum FWHM « 3000 + 1000 k m s 1 
(Schmid et al J 120011) . The optical/UV flux is highly polarised 
(Martin = et = al. 1983). VLT-spectropolarimetry measurements by 
ISchmid et al. ( 200 1) revealed a polarisation degree ranging from 
5% (red) to 13% (UV), which are one of the highest polarisation 
degrees observed for type 1 objects. 

Fairall 51 is very bright in X-rays (about 1 millicrab). It was 
observed twice with high-quality X-ray spectrometers on-board 
the XMM-Newton satellite. The first observation was performed 
in 2005 with the observed flux Tjjfg kcV ~ 3 x 10 _11 erg cm -2 s -1 , 
and the second one about half a year later with the observed flux 
■^ 2 -TokeV ~ 0-9x 10~ n erg cm~ 2 3 s _1 . The spectra significantly dif¬ 
fer and the complex soft X-ray spectrum may be described by the 
presence of three zones of a partially-covering warm absorber 
with high, moderate and very low ionisation dRicci et HIM. 

Short-time variability of this source was suggested from 
two Swift observations performed in 2008. The observations 
were separated only by five days. Although the exposures were 
very short (a few ks each) the spectra were qualitatively dif - 
ferent below 7 keV dRicci et alJ 120101 : IBeuchert et alJ 120131) . 
iBeuchert et alJ (1201 3h determined that the variable absorber must 
be located not further than the Broad Line Region based on the 
short variability time-scale. 

In this paper, we present the results of a recent monitoring of 
Fairall 51 by the Japanese X-ray satellite Suzaku dMitsuda et al.l 
l2007h , during which four 30 ks exposures were performed sepa¬ 


rated by 1.5, 2 and 5.5 days. The paper is organised as follows: 
the observations and the data reduction are described in Sect. [2] 
Results are presented in Sect. [3 and discussed in Sect. [4] Main 
conclusions are summarised in Sect. 0 

2. Observations and data reduction 

We observed Fairall 51 using the X-ray Suzaku satellite 
dMitsuda et al.ll20Q7h during four ~30 ks long exposures in the 
first half of September 2013 (see Table[T|). The observation iden¬ 
tification numbers are 708046010, 708046020, 708046030, and 
708046040, respectively. All the observations were performed at 
the XlS-nominal pointing position. 

The Heasoft package version 6. l43 was used for the data re¬ 
duction and also for the subsequent spectral and timing analysis. 
The data were processed standardly following the Suzaku Data 
Reduction Guid^l (version 4). For all XIS detectors, we com¬ 
bined both 3x3 and 5x5 modes to extract the event files. The 
source spectra were obtained from a circle around the centre of 
the point spread function with the radius of 260". We defined the 
background extraction region as an annulus around the source 
circle with the outer radius of 360" to avoid any contamination 
from the calibration source at the edges. We created related re¬ 
sponse matrices and ancillary response files using the tools xis- 
rmfgen and xissimarfgen. The HXD/PIN spectra were reduced 
with the tool hxdpinxbpi. The tuned background files0 were used 
to model the non X-ray background. The cosmic X-ray back¬ 
ground was estimated using the model bv lBoldtl (1987) in accor¬ 
dance with the Suzaku Data Reduction Guide. 

We used the Xspec software (fArnaudll 19961) version 12.8.2 
for the spectral analysis. The cross-normalisation factors be¬ 
tween the spectra of different Suzaku detectors were free for 
XIS 1 and XIS 3, and fixed to 1.16 for HXD/PIN according to the 
Suzaku Data Reduction Guide. Fits were performed in the 0.5- 
lOkeV energy interval for the XIS spectrum, and from 15keV 
for the HXD/PIN spectrum. We fit the HXD/PIN spectra up to 
the energy corresponding to a 5% source detection level above 
the background level which is well above the reported systematic 
uncertainty in the Suzaku Data Reduction Guide. All HXD/PIN 
spectra but the second one were considered in the 15-55 keV 
energy range. Only the HXD/PIN spectrum of the second ob¬ 
servation was limited by 45 keV due to a lower signal-to-noise 
ratio caused by a shorter exposure time («22ks) compared to 
the other observations. 

We used C-statistics ( Cashlll979l) for fitting the data. How¬ 
ever, we also express the goodness of the fit with more familiar 
X 2 values that were used as a test statistics. Only for this purpose, 
we binned the spectra to contain at least 30 counts per bin. All 
the quoted errors correspond to a 90% confidence level for one 
interesting parameter. 

3. Results 

3.1. Timing properties 

We performed a series of four observations separated by 1.5, 2 
and 5.5 day intervals to constrain the characteristic time scale 
of the short-term variability of Fairall 51. Figure [T] (left panel) 
shows the XIS0 and HXD/PIN light curves from all four obser¬ 
vations. While the soft X-ray flux decreased by a factor of two 

1 http://heasarc.nasa.gov/lheasoft/ 

2 http://heasarc.nasa.gov/docs/suzaku/analysis/abc/ 

3 We used a new version 2.2 updated on 2014 Jun 9 
(http://heasarc.gsfc.nasa.gov/docs/suzaku/aehp_data_analysis.html). 
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Table 1. List of Suzaku observations from September 2013. 


Observation 

Identification Number 

Start Time 

Stop Time 

Net Exposure [ks] 

1 

708046010 

Sep 4th 01:45:56 

Sep 4th 12:31:08 

31.5 

2 

708046020 

Sep 5th 19:22:38 

Sep 6th 03:59:00 

31.0 

3 

708046030 

Sep 7th 19:35:09 

Sep 8th 10:17:13 

24.4 

4 

708046040 

Sep 13th 12:16:30 

Sep 14th 07:32:52 

30.4 
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Fig. 1. Left: XIS 0 (upper, with small circles) and HXD/PIN (lower) (all background subtracted) light curves in the 0.5 - lOkeV, or 15 - 55 keV, 
energy range, respectively. Right: The XIS 0 light curves in two different energy bands 0.5-3 keV and 3-10 keV (all background subtracted), and 
the hardness ratio between the hard and the soft band. The gaps between the observations are shrunk for the clarity. The time bin size is 2880 s. 


between the first and the last observation, the hard X-ray flux 
stays more or less constant. 

The soft X-ray variability is shown in more details in the 
right panel of Fig. [T| where the light curves in two energy bands, 
0.5-3 and 3-10 keV, are shown together with the hardness ratio 
defined as the flux ratio between these two bands (hard to soft). 
The major change in the hardness is measured during the last ob¬ 
servation, when the hardness ratio is significantly larger than in 
the previous three observations suggesting a significant change 
of the spectrum during this observation. 

The flux significantly varied also within the individual ob¬ 
servations. It decreased almost to its half during the second ob¬ 
servation, but the source became brighter again during the third 
observation. This variability is accompanied with a small hump 
in the hardness ratio. However, a low signal-to-noise ratio limit 
prevented us from a detailed investigation of this feature in the 
light curve, and we further deal with four time-averaged spectra 
of the individual observations. 

3.2. First took at X-ray spectra 

To have a first look at the spectra, we employed a simple model 
consisting of a power law with its normalisation as the only vari¬ 
able parameter between the observations. The power law was 
affected by the photoelectric absorption by cold material, for 
which we used the TBABS model (iWilms et alil200d ) with So¬ 
lar abundances. The column density was a free parameter but 
linked between the four observations. 

Figure [2] shows the X-ray spectra of all Suzaku observations 
together with the data residuals from this simple model. The 


power-law photon index is T ~ 1.7 and the absorber’s column 
density is Ah ~ 3.8 x 10 22 cm 2 . The quality of the fit is not 
acceptable with C/v = 12033/6454 « 1.9 and with clear data 
residuals from the model over the entire energy range (see the 
left bottom panel of Fig. [2}. 

The most evident spectral variability is at 2-6 ke V, while the 
spectra do not seem to vary at very high and also very soft (< 
1 keV) X-ray energies. The latter suggests that the spectrum is 
dominated by two different spectral components in the 0.5-1 and 
2-6 keV energy bands. A dominating spectral component at the 
very low energies can be a scattered power-law emission that 
is not affected by a variable absorption of a circumnuclear mat¬ 
ter that strongly affects the direct emission, similarly to what 
has been found for other similar sources, e.g., ESO 323-G77 
dMiniutti et al.ll2014l) . 

Other data residuals from the simple power-law model occur 
at the iron line band (see the right panel of Fig. [2}. The emission 
line profile itself seems to be asymmetrically broadened towards 
its red wing, especially in the lowest-flux observation. However, 
it would be premature to draw any conclusions about the iron 
line before finding an appropriate description for the continuum. 
Therefore, we start with a Gaussian model for the iron line. Be¬ 
yond the emission iron line, two absorption features are visible 
in the spectra. 

The absorption features above the iron line are likely due 
to an ionised gas in the line-of-sight, probabl y ass ociated to the 
warm absorber that was already revealed by iRicci et akl (1201 Of) 
who identified in total three ionised absorbers with different ion¬ 
isation based on the XMM-Newton spectral analysis. 
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Fig. 2. Left: Up: X-ray spectra of four Suzaku observations (1st black, 2nd red, 3rd green, 4th blue). Bottom: Data residuals from a simple 
absorbed power-law model with T ~ 1.7. Right: A more detailed look at the iron line and Compton hump as residuals from the simple absorbed 
power-law model. Only XIS 0 and HXD/PIN spectra are plotted for clarity. Data are re-binned for plotting purposes only. 
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Fig. 3. Data residuals from the different models for all observations (1st 
black, 2nd red, 3rd green, 4th blue). From top to bottom: Residuals 
from the model A-D. 
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The simple power-law model fails to describe the hard X-ray 
emission (see the right panel of Fig. |2}. The spectral curvature 
might be explained by a reflection hump. Its presence close to 
the lower energy limit of the HXD/PIN spectrum suggests that 
the reflection emission is red-shifted. 

Finally, the spectral features around 2 keV might be affected 
by calibration uncertainties and therefore we further omit the 
data at 1.6-2.4keV. 

Summarising the results of the first-look analysis, Fairall 51 
is variable mainly at the 2-6 keV energy band while there is only 
little variability in the hard X-rays and also around 1 keV. 

3.3. Modelling the complex absorption 

Next, we employed a model with two power-law components 
with the photon index linked between them - one affected by the 
local absorber and one not. The local absorber was first modelled 


by a cold absorber using ZTBABS model with the column den¬ 
sity free to vary between the observations. We further used two 
Gaussian lines to model the absorption features in the iron-line 
energy band, with free energies and the redshift fixed to the cos¬ 
mological redshift of the galaxy. Furthermore, we used a Gaus¬ 
sian line profile to account for the emission iron Ka line. The 
global model was then absorbed by the interstellar absorption 
in our Galaxy with the column density A'n = 6.9 x 10 2l) cm~ 2 
dKalberla et alJ[2005) using the TBABS model. Solar abundances 
were assumed in all absorption models. 

We obtained a statistically better fit but not yet acceptable 
with C/v = 7730/5798 w 1.3. The column density of the vari¬ 
able absorber was found to increase from Nn = 2.95/Jj ^ x 10 22 
cm -2 during the 1st observation to Nu = 6.8 ± 0.3 x 10 22 cm~ 2 
during the 4th observation. For the iron K a emission line, we 
obtained E = 6.32 + 0.03 keV, <x = 0.14 ± 0.03 keV and equiv¬ 
alent width ranging from EW = 120/™ eV (1st observation) to 
EW = 320 ± 40 eV (4th observation). The line was found to be 
broadened and red-shifted from its intrinsic 6.4keV energy. The 
FWHM widths of the absorption lines were found even larger. 
Their values were pegged at the upper limits, which we set to 
cr = 0.3 keV. The energy of the first absorption line was found 
to be £ = 6.75 ±0.1 keV for the first three observations and 
E = 6.6 ± O.lkeV for the 4th observation. The energy of the 
second absorption line was found to be £ = 8.2 ± 0.2 keV for all 
observations, consistent within the uncertainties. 

The power-law photon index in this simple model was 
r ~ 1.5. Unusually flat spectral slopes have been often found 
to be due to a part i ally-covering ionised absorber (see e.g . 
iPiconcelli et al.ll2004t iMathur et al .112009k ISvoboda et al.ll2012ft . 


As a next step, we therefore employed a partially-covering 
ionised absorption to model the absorption features in the spec¬ 
tra. We replaced the absorption Gaussian lines by the ZXIPCF 


model (Reeves et al. 

2008 

) using the XSTAR photo-ionisation 

code (Kallman & Bautista 

2001). We first start with one com- 


ponent and refer to this model as “model A”. The best-fit values 
are presented in Table [2] and the data residuals from the model 
are shown in the upper panel of Fig. [3] The photon index of the 
power law F » 1.8 is consist ent with the characteristic slope of 
Seyfe rt 1 galaxies (see e.g. iBianchi et al.ll2009at iMalizia et~aTl 
l2014i) . We found high ionisation of the absorber, log£ > 3, con¬ 
sistently with iRicci et aD ( 2010 ). but with the column density 
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Table 2. Spectral analysis of Fairall 51 - model parameters. 


Model Component 

Parameter 

Model A 

Model B 

Model C 

Model D 



1st obs 

4th obs 

1st obs 

4th obs 

1st obs 

4th obs 

1st obs 

4th obs 

power law 

r 

1 84+ 003 

1,0 -0.02 

1.80 ±0.02 

1 96+ 004 

A -^ u -o.oi 

2.26 ± 0.04 


X direct [10- 3 ] 

11.5/?;/ 

10 6 +1 ° 
au ' u -0.4 

11 o +0 - 3 

1 a,u -0.5 

Q 0+0.5 
—0.4 

13 7 +10 

AJ *'-0.3 

11 o +0 - 8 

A a,u -0.3 

20 7 +L4 
^ u -' -2.0 

10 i +0 - 7 

au,a -0.9 


^scattered [10 ] 

0.19/°;°/ 0.17 ±0.01 

0.18/°;°/ 0.17 ±0.01 

0 12 +002 
u - iZ, -0.03 

0 13 +002 
u - A:, -o.oi 

0.13 ±0.02 0.14/°;°/ 

local cold absorber 

Nh[ 10 22 cm- 2 ] 

3.9 ±0.1 7.6 ±0.4 

o 7+0.2 
-0.1 

7 1 +0.4 
'-0.3 

- 

- 

low-ionised absorber 

N h [ 10 22 cm- 2 ] 

- 

- 

2-7!°: 5 0 5.7 ± 0.4 

7 q+0.7 
—0.4 

6 7 +u 


log ft 

- 

- 

1 l +0 1 

1 ’ 1 -0.2 

1.2 ±0.1 

mid-ionised absorber 2 

N h [ 10 22 cm- 2 ] 

- 

- 

7 3 +L3 

' -0.5 

8* 2/4*2 

7.7 ± 0.4 6.5/J f 


log ft 

- 

- 

1.9 ±0.1 

1.6 ±0.1 

high-ionised absorber 

N h [ 10 22 cm- 2 ] 

26/® 114 ±17 

18!’ 

98/5 

23 +15 

^ J -10 

28// 4 


log ft 

3.3 ±0.1 

a 7+0.1 
d ‘-0.2 

3 6 +0 - 2 
• z, - u -o.i 

3.6 ±0.2 


Voutdow [kms -1 ] 

3300 ± 600 

4200 +130 ° 

t +z,L/u_i4oo 

0 ± 1500 

1400//°° 

very-high ionised absorber 

A h [ 10 22 cm- 2 ] 

- 

23+ 27 
^ -10 

8// 3 

23/3 

- 


log ft 

- 

4.2 ±0.2 3.6-:;' 

A oq+0.09 
-0.15 

- 


^outflow [km s ] 

- 

300/400° 

27QO+ 800 
^ u -1600 

- 

Gaussian line 

£[keV] 

6.37 ± 0.02 

6.38/°;° 3 

- 

- 


<r[keV] 

0 03 +0 04 

-0.03 

0 04 +0 06 

u - w ^-0.04 

- 

- 


normalisation [ 10 ^ ] 

1.2 ±0.2 

1 2 +05 

1 ‘ -0.3 

- 

- 


equivalent width [eV] 

35// 

561‘‘ 

36// 6 

58!? 4 

- 

- 

relativistic reflection 

spin a/M 

- 

- 

- 

0.80/°;!/ 


inclination [deg] 

- 

- 

- 

20 +6 
^ -20 


radial emissivity 

- 

- 

- 

3.5 ±0.3 


normalisation [lCL 3 ] 

- 

- 

- 

54//° 

71+25 
' -30 

distant reflection 

iron abundance [Solar] 

- 

- 

< 0.35 

0 89+ 018 
u -°^-0.14 


normalisation [ 10 ’ ] 

- 

- 

6 3 +1 ° 
u -0.6 

Q 4+2.5 

fit goodness 

C/v 

6399/5809 

6297/5801 

6094/5810 

5988/5806 

test statistics 

X 2 /v 

6410/5809 

6294/5801 

6097/5810 

5967/5806 


Notes. The values of the model parameters which were allowed to vary between the four observations are reported for the 1st (highest-flux) and 4th (lowest-flux) observations only for brevity. A 
single value for all observations is reported when the parameters were free in fitting but fixed to the same value or did not significantly change between the observations. All the errors correspond to 
the 90% confidence level for one interesting parameter. All normalisations are expressed in units of photons keV~* cnr 2 s -1 at 1 keV. 
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Fig. 5. The final model (model D) in the individual observations (1st 
black, 2nd red, 3rd green, 4th blue). 

A'h ~ 10 23 - 10 24 crrr 2 , roughly an order of magnitude larger. 
However, the global fit is not satisfactory with C/v ~ 1.1, and, 
in particular, residuals around ~ 7 keV indicate the need for an¬ 
other high-ionised absorption component. 

The addition of another warm absorber component signifi¬ 
cantly improved the fit (from C/v = 6399/5809 to 6297/5801). 
The statistical F-test using the y 1 test statistics gives the proba¬ 
bility « 3 x 10 19 that the improvement is just coincidental. We 
refer to this model as “model B” and the best-fit values are sum¬ 
marised in the Table [2] The value of the ionisation of the second 
warm absorber is higher, log £ w 4, and it accounts for the resid¬ 
uals above 7keV (see the second panel of Fig. [3j. We allowed 
also the covering factors to vary, and we tried several tests with 
different initial values, but we found no statistical improvement 
compared to fits when these parameters were fixed to 1. The ion¬ 
isation parameter was also left free but no significant change was 
detected across the observations. Only column densities of both 
cold and warm absorber were required to vary significantly. 

As a next step, we allowed the cold absorber to be ionised 
and partially covering. In XSPEC, we replaced the ZTBABS 
model by another ZXIPCF component. We first allowed also the 
redshift of this absorber to vary. However, its value reached 0.1, 
which would correspond to an inflow of > 25000 km s _1 . This is 
a quite extreme value that has not been reported in any similar 


source. Because the redshift in this case is not obtained from the 
frequency shift of a clearly defined absorption line, we consid¬ 
ered this measurement as not indicative and we rather fix the red¬ 
shift to the cosmological value of the galaxy. The obtained fit is 
significantly better than the previous one with C/v = 6228/5804, 
i.e., A C/v = 69. However, we realised that instead of having two 
highly-ionised and one low-ionised component, the fitting pro¬ 
cedure converged to a low-ionised, mid-ionised and one highly- 
ionised component. Contrary to the previous model, the column 
densities of the high-ionised warm absorber were not required 
to vary between observations. Also, we found that the ionisation 
was not required to vary, and covering factors were consistent 
with being 1 for all observations. As a check of this result, we 
linked the column densities, and checked whether variable cov¬ 
ering factors and ionisations may describe the data as well. How¬ 
ever, we obtained statistically worse fit with C/v = 6259/5795. 

We checked the residuals and realised that the most signif¬ 
icant improvement compared to model B appeared at energies 
< 1 keV, while the residuals at the iron line energy band got 
worse. We therefore added another absorber to account again 
for high ionisation. We found a new statistically better fit with 
four absorbers (C/v = 6138/5797, A C/v = 90 from the previous 
model). However, we did not obtain any significant improvement 
in modelling the residuals at energies above lOkeV. 

3.4. Modelling the reflection 

The models presented so far were not able to fully describe 
the data residuals in the iron line band and around 15keV. 
The presence of iron emission line is required in all mod¬ 
els, which indicates a significant contribution of X-ray reflec¬ 
tion to the observed spectrum. To consistently treat the X- 
ray reflection, we replaced the G auss ian line by the reflection 
model PEXMON (iNandra et al.1 l 2007h that combines the iron 
line and the reflection continuum using the PEXRAV model 
(iMagdziarz & Zdziarskil 19951) . We used the PEXMON model 
only for reflection, i.e., with the reflection fraction parameter R 
fixed to -10 The inclination was fixed to 45 degrees. We refer 
to this model as “model C” and the best-fit values are presented 
in the Table [2] The statistical improvement compared to the pre¬ 
vious model is A C/v = 43. However, the addition of the distant 
reflection component still did not improve satisfactorily the most 
prominent residuals at the HXD/PIN energy band (see the third 
panel of Fig.|3}. 

As a next step, we employed reflection from the innermost 
region of the accretion disc where the radiation is significantly 
red-shifted due to effects of the general relativity. We therefore 
added a second reflection component to the model that is con¬ 
volved with a relativistic disc kernel dDovciak et al .1120041) . Free 
parameters of the relativistic model were the black hole spin, in¬ 
clination angle and radial-emissivity index. The inner radius was 
set to the marginally stable orbit depending on the spin value. 
The outer radius was set to 400 r g , where r g = Of is the gravita¬ 
tional radius. The iron abundances were linked to the abundances 
of the distant reflection model. We still kept the inclination of the 
distant reflection fixed to 45 degrees, because the global fit is not 
very sensitive to this parameter. 

We further assumed that the relativistic reflection is affected 
by all absorbers while the neutral reflection coming possibly 

4 Only to estimate the reflection strength we also used the PEXMON 
model for primary power-law emission. We fixed the power-law nor¬ 
malisation to zero and let the parameter R to vary. We got the reflection 
strength R ~ 0.5. 
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Fig. 6 . Relation between the intrinsic 2-10 keV flux and the column 
density of the cold (least ionised) absorber for the model D (black, cir¬ 
cles) and model B (red, triangles). 


from the distant polar scattering region is unaffected by the 
warm absorber. We obtained a significantly improved fit (C/v = 
5988/5806, A C/v = 106 from the previous model). The best- 
fit values are presented in the Table[2]as model “D” parameters. 
We considered only three absorbers in this model, because we 
realised during the fitting procedure that the highest ionised ab¬ 
sorber becomes less important with the new model (the differ¬ 
ence between the model with three and four absorbers was only 
AC « 10). Some of the absorption features at 7-10 keV were 
likely attributed to absorption edges of the reflection models. 

The final model is statistically well acceptable with no sig¬ 
nificant residuals across the entire Suzaku bandpass (see the bot¬ 
tom panel of Fig. [3}. The model with its components is shown 
in Fig. [4] Figure 0 shows the comparison of the final model in 
all four observations. The low-ionised absorber is the most re¬ 
sponsible for the curvature of the spectrum at energies 1-10 keV. 
The scattered power-law emission dominates at 1 keV over the 
absorbed intrinsic flux and implies that the flux is not strongly 
variable there. The hard energy band is shaped by reflection. Rel¬ 
ativistic smearing is responsible for shifting the Compton hump 
to lower energies. 

We note a minor discrepancy between the data residuals and 
model predictions in the lower edge of the XIS spectrum, com¬ 
pare Fig. [2] and Fig. 0 at 0.5-0.7keV. This can be, however, an 
instrumental effect due to the redistribution of the more ener¬ 
getic photons to lower energies that appears at the low edge of 
the spectrum rtKovama et aljl2007h . 


4. Discussion 

4. 1. Variable spectral components 

Our monitoring observational programme revealed that the 
short-term variability of Fairall 51 is partly caused by the vari¬ 
able column density of the absorbing gas, and partly by the 
changes in the intrinsic luminosity that also significantly varies 
in the soft X-rays between the observations. Table [3] shows our 
measured values of the observed X-ray flux in the 2-10 keV 
range, the intrinsic 2-10 keV fluxes constrained from our final 
model, the observed luminosity determined using the cosmolog¬ 



Fig. 7. Contour plots between the column densities of the variable low- 
and mid-ionised absorbers. The colours correspond to different observa¬ 
tions (1st black, 2nd red, 3rd green, 4th blue). The contours correspond 
to lcr, 2<r, and 3cr levels. The best-fit values are marked by a small cross 
at the graph. 

ical distanc^B, and the estimation of the bolometric and ionising 
luminosity. The bolometric luminosity was_constrained from the 
X-ray luminosity using a relation bv lMarconi et al.l d2004l) . 

Because the intrinsic X-ray luminosity decreased almost by 
a factor of two according to the final model, we get different val¬ 
ues for the bolometric luminosity. We note, however, that these 
are only estimates based on measured X-ray luminosities, and 
the bolometric luminosity itself does not need to vary. The es¬ 
timated values were Lb 0 i ~ 3.7 x 10 44 ergs _l for the first ob¬ 
servation, and Lboi ~ 2 x 10 44 erg s _1 for the last observation. 
This translates to Eddington ratios Tedd ~ 0.03, and /fedd ~ 0.02, 
respectively. We used relation for the Eddington luminosity as 
^Edd = 1.26 x 10(A7/;V/ 0 ) e rg s 1 and the value of the mass 
M » 10 8 M 0 dBennert et alJfeoOfil) . 

To constrain the ionising luminosity, we used the rela¬ 
tion by lYasudevan et al.l d2010l) to estimate the soft luminosity 
Ti 3 . 6 -iooeV = ^13.6-iooevTboi- where /ci3.6-iooeV goes from 0.21 
to 0.59 for accretion rates from 0.01 to 0.61. By the linear in¬ 
terpolation, we constrained the factors for each observation de¬ 
pending on the estimated accretion rate. We determined the X- 
ray luminosity L 0 . 1 - 13 . 6 kev by an extrapolation of the measured 
X-ray luminosity in 0.5-10keV. To get the ionising luminosity, 
we summed the soft and the X-ray luminosity. The results are 
shown in the last column of Table[3] Figure [6] shows the relation 
between the column density of the least-ionised absorber and the 
intrinsic 2-10 keV flux. This relation depends on the employed 
model, and therefore we show this relation for models B and D. 
Both results strongly indicate that the column density of the ab¬ 
sorber increases with decreasing intrinsic flux. 

Surprisingly, we have not detected any significant variabil¬ 
ity in the ionisation of the absorbers that would reflect the de¬ 
crease of the illuminating flux. This would suggest that the ab¬ 
sorber is not in the photo-ionisation equilibrium. Such behaviour 
has been reported also for some other so urces, such as Mrk 33 5 
dLonginottiet al.ll2013l) . or NGC 3783 (iKrongold et al.l l2005h . 
We note, however, that in our case, the models of the low- and 
mid-ionised absorbers are based on the continuum spectral shape 
and not on the individual absorption lines that would define 

5 using the cosmological parameters H 0 = 70km s -1 Mpc~*, q 0 = 0, 
and do = 0.73 
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Table 3. Observed and intrinsic flux and luminosity measurements of Fairall 51. 


observation 

number 

f*2-10 keV,observed 

[10 11 erg cm -2 s~*] 

b2-10keV, intrinsic 

[10~ u ergcnr 2 s -1 ] 

^2-10 keV,observed 

[ 10 43 erg s -1 ] 

L1 5 -50 ke V,observed 

[10 43 erg s _1 ] 

Lbolometric 

[10 44 erg s _1 ] 

L'O.0136-13.6 keV,intrinsic 

[10 44 erg s 4 ] 

i 

2.3 ±0.1 

4.1 ±0.4 

1.8 ±0.2 

1.9 ±0.2 

3.7 ± 0.4 

1.5 ±0.2 

2 

2.1 ±0.1 

3.6 ± 0.4 

1.6 ±0.2 

1.9 ±0.2 

3.2 ± 0.4 

1.3 ±0.2 

3 

1.8 ±0.1 

3.1 ±0.3 

1.4 ±0.1 

1.8 ±0.2 

2.7 ± 0.3 

1.1 ±0.1 

4 

1.3 ±0.1 

2.5 ± 0.3 

1.1 ±0.1 

1.7 ±0.2 

2.0 ±0.3 

0.8 ±0.1 




Fig. 8. Contour plots between the column densities of the variable cold and warm absorbers. The colours correspond to different observations (1st 
black solid, 2nd red dashed, 3rd green dotted, 4th blue dash-dotted). The contours correspond to lcr, 2<x, and 3cr levels. Left: Results obtained 
with the model B. Right: Results obtained with the model D with highly-ionised warm absorber allowed to vary. 


the ionisation parameter more precisely. Therefore, no measured 
variability of the ionisation parameters may be simply due to 
coarse resolution of the CCD detectors. 

All the employed models suggest that the most variable spec¬ 
tral component is the least-ionised absorber. Figure |7] shows a 
contour plot between the column densities of the low- and mid- 
ionised absorbers in the final model D. All the free parameters of 
the final model were allowed to vary during the contour calcula¬ 
tions. The column density of the low-ionised absorber changed 
by a factor of two between the first three and the last observation, 
while the column density of the mid-ionised absorber may stay 
unchanged within 3<x levels. 

As the values of column densities are similar and the ion¬ 
isation parameters are not constrained from any discrete fea¬ 
tures, we tested a scenario, in which the low- and mid-ionised 
absorbers are replaced by one variable absorber. We considered 
two cases: (1) the column density and the ionisation are vari¬ 
able between the observations, and (2) the column density and 
the covering factor are variable. In both cases, the resulting fit 
is worse by AC » 30 compared to the final model. The F-test 
using x 1 statistics gives the probability pa 3 x 1CL 6 that the 
statistical improvement with two absorbers is just coincidental. 
Higher-quality data would be needed to more precisely constrain 
the structure of the absorbing gas. 

Further, we investigated possible variability of the highest 
ionised absorber. Figure [8] shows contour plots between the col¬ 
umn density of the cold and highly ionised absorber. The left 
panel shows results obtained by the model B. The right panel 
shows results obtained by a modified version of the model D, 
in which we allowed the warm absorber with high ionisation 


to vary instead of the mid-ionised absorber. While the warm- 
absorber column density is required to change in the model B, 
the model D shows that this parameter may stay unchanged dur¬ 
ing all observations. 


4.2. Constraints of the variable absorber’s location 


We found that the spectral change between the last observation 
from the previous three observations is mainly due to the in¬ 
crease of the column density of the cold or least ionised absorber. 
Based on the measured time scale of the spectral variability, we 
may estimate the location of this variable absorber. Let us as¬ 
sume that the main velocity of the absorbing cloud is the orbital 
velocity around the centre that can be approximated as the Kep- 
lerian velocity at the distance R: 


vk 


G/V/bh 


R 


( 1 ) 


where Mbh is the black-hole mass, and G = 6.67 x 

10 -8 cm 3 g -1 s~ 2 is the gravitational constant. For a cloud with 
characteristic size s, we can write: 

At 2 

R = GM bh —, (2) 

s l 

where At is the time difference between the start and the end of 
occultation. Assuming a spherical symmetry and homogeneity, 
we can estimate the size of the cloud from the measured column 
density Ah: 


S a 


Ah 

n 


(3) 
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where n is the particle density. We can now rewrite the eq.[2]as: 
A t 2 n 2 


R « GMbh 


ANi 


(4) 


where ANh is the change of the column density in the time inter¬ 
val At. We may rewrite the eq.[4]to have a form more appropriate 
to the typical values that are measured. We obtain an equation 
analogous to the eq. 3 in iRisaliti et al.1 (12002!) : 


R * 1.335 x 10 


13 


M B h 

10 7 M f 


f H f 

\10 6 cm” 3 / 


AN h 


At 


1 Ms / \ 10 22 cm” 2 


(5) 


Using this equation, we can constrain the distance of the ab¬ 
sorber from the change of the column density in a given time 
if we know the mass and the volume density. The mass of the 
black hole is usually constrained from different methods in¬ 
cluding reverbe ration mapping (Petersonl 120141) . stellar veloc- 
ity dispersions (iF errarese & Merritt]' 2000b, bulge lu minosities 
(iKormendv & Richstonel 1995T Magorrian et al.l 1 19981) . or from 
single optical emission lines due to a ti ght correlation between 
the luminosity and mass (see, e.g.. IShen et al.il20(M and refer¬ 
ences therein). However, the volume density is unknown. The 
density should not exceed 10 6 particles per cm ’ if it is part of a 
torus and it should be about 10 9 particles_per cm 3 if the cloud 
belongs to the Broad Line Region (iWeedmanll 19771) . 

The volume density can be estimated if the ionisation £ of 
the absorber is measured: 


L 

fR 2 ’ 


( 6 ) 


where L is the illuminating luminosity the energy range 13.6 eV 
- 13.6keV. We can now replace the density in eq.|4]and get: 


R 


GMbh 


AfL 2 

AAfrF 


(7) 


Again, we can rewrite this equation with the quantities nor¬ 
malised to their typical values and we get an equation analogous 
to the eq. 3 in lLamer et al.1 d2003l) : 


R a 2.66 x 10 17 
L 


M B h 

<N 

< 

AN h \ 

10 7 M 0 

[lMs) 

10 22 cm -2 / 


-2 


10 43 erg s~ 


erg cm” 2 s” 1 


- 2 . 


( 8 ) 


The black-hole mas s of Fairall 51 was estimated as M » 
10 8 M g by iBennert et alJ (120061) based on the relation between 
the luminosity at 5100 A and the black hole mass (i Peterson et al.l 
l2004i) . We measured the change of the column density of the 
cold absorber ANh ~ 4 x 10 22 cm” 2 to be * 0.45 Ms (see Fig.Q}. 
The ionising luminosity is taken as an average value between the 
3rd and 4th observation, L ~ 10 44 ergs”’ (see Table [3}, and the 
ionisation is £ « 15 erg cm” 2 s” 1 (see Table [2}- After inserting 
these values into the eq.[8]we get the estimation for the distance 
of the absorbing cloud from the centre: 

R ss 1.5 x 10 l7 cm ~ 0.05 pc ~ 60 light days. 


The distance of the Broad Line Region can be estimated from 
the optical luminosity. A tight relation between the BLR radius 


and th e luminosity at 5100 A w as constraine d by Kaspi et al. 
(120051) and recently updated by iBentz et akl (120131) who also 
subtracted contributions of the host galaxies. The luminosity of 
Fairall 51 is L 5 , (I()A ~ 10 43 erg s” 1 iBennert et al.l 12006 ). This 
corresponds to the BLR size of ~ 40 light days according to 
iBentz et akl (120131) . This value is consistent what we found for 
the variable X-ray absorber, taking into account the general un¬ 
certainties of both estimations. 

We can also derive the density of the absorbing cloud from 
our measurements. Using eq. [6j we obtain: 

n w 3 x 10 s cm” 3 . 

This value is slightly lower than the averaged expected density 
of BLR clouds dPetersonll2014l) . Our measurements of slightly 
lower density, larger radius and low ionisation are consistent 
with the reverberation in optical spectral lines that revealed a 
stratification in the ionisation and the density of the BLR region 
(see, e.g. JGaskell & Snarkeil 19861) . both decreasing with the ra¬ 
dius. 

Although our measured location of the variable absorber is 
fully consistent with the Broad Line Region, we note that the in¬ 
nermost part of the torus cannot be ruled out either, as its value is 
not well established. The most recent mid-infrared observations 
of nearby galaxies revealed that there is a large scatter in the ge¬ 
ometrical properties of their dusty tori (iBurtscher et al.l 120131) . 
and even very complicated structures of the dust within the cen¬ 
tral parsec were reported dTristram et al.ll2014l) . 

4.3. Scattering region 

The presence of a scattering region that is not located on the 
line-of-sight is suggested from the optical polarisation measure¬ 
ments. Our X-ray spectra show significant spectral variability in 
the energy range 2-6 keV (see Fig. 0 that is best explained by 
a variable absorption. Less variability is evident around 1 keV, 
which can be well explained by a scattered power-law emission 
that is not affected by a variable absorber. Although its intrin¬ 
sic normalisation is lower than the normalisation of the direct 
power-law emission by two orders of magnitude, it dominates 
the spectrum at * 1 keV where the direct nuclear emission is 
heavily absorbed (see Fig. |4]). 

The scattering region could be associated with the extended 
Narrow Line Region located at a rotational axis with the sim¬ 
ilar geometry as considered by iMiniutti et al.l d2Q14t) for an¬ 
other polar-scattered Seyfert 1 galaxy ESO 323-G77 (see their 
Figure 8), or in obscured Seyfert type 2 galaxies (see, e.g., 
iBianchi et alj|2006l) . This region is also supposed to be respon¬ 
sible for the measured polarization degree that is unusually high 
for a type 1 object dSchmid et al.l200lf) . 

An alternative explanation for the high polarisation degree 
in polar-scattered Seyfert 1 galaxies might be due to scattering 
on disc winds (on the external layer s of the absorbing clouds) a t 
the opposite s ide from the observer dMarin & Go osmann 20131) . 
|m arinl (120141) showed on an AGN sample including Fairall 51 
that scattering on the disc wind better corresponds to the data 
than scattering in a polar region. Nevertheless, to draw a clear 
conclusion about the geometry of the scattering region. X-ray 
polarimetric measurements would be required. 

We note that two components were also detected for the re¬ 
flection. A part of reflected emission is coming from the inner¬ 
most accretion disc and is also affected by the absorber. The 
other part originates at a farther region and is not smeared by 
the relativistic effects. The distant reflection may occur at the 
same scattering region as for the power law. 
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Fig. 9. Contour plots between the spin and the inclination in the rela¬ 
tivistic model. The contours show lcr, 2<x, and 3<x levels. 


4.4. Geometrical constraints from the reflection models 


The reflection features have been detected mainly owing to the 
broad-band coverage of the Suzaku satellite. Figure [2] shows the 
presence of both the iron line and the Compton hump. Although 
the very simple power-law model revealed a skewed profile of 
the iron line especially for the 4th observation (see the right 
panel of Fig. 0, these data residuals could be reasonably well 
modelled by a complex absorption and distant reflection (see 

Fig-0- 

Flowever, the distant reflection was not able to reproduce 
the data residuals at energies around 15keV. We found that the 
Compton hump is significantly red-shifted, which can be well 
accounted by relativistic reflection. Although hard excesses have 
been successfully expl ained also by other models, incl uding ab¬ 
sorption scenario (e.g. jMiller et al.||2008 ; Tatum et all 201 3l) or 
complex reprocessing in the torus (Murphv & Yaqoob 2009), we 
have not succeeded with simple tests using these alternatives to 
describe the data. We leave a more complex analysis using these 
alternatives for future work. 

Using the relativistic reflection model, we were able to mea¬ 
sure the spin and the inclination angle of the accretion disc, 
though with a large uncertainty given the complexity of the fi¬ 
nal model. We obtained a ~ 0.8, and i < 20°. The contour plot 
between these two parameters is shown in Fig.[9](the model pa¬ 
rameters that are shown in Table [2] were all left free during this 
contour calculations). The spin is not constrained sharply. Its 
value can be any between 0.6 and the maximum value within 
the 3cr level. The inclination is lower than it would be expected 
from the polarisation measurements that suggest the inclination 
to be around 45 degrees. A possible explanation could be that 
the relativistic reflection also comprises a direct and scattered 
component. The scattered part can be relatively enhanced if the 
direct one is heavily obscured. However, in the final model, the 
scattered component of the primary power law is two orders of 
magnitudes lower than the direct part. Assuming that the primary 
and disc reflection emission do not originate far from each other, 
the explanation by reflection dominated by a scattered compo¬ 
nent seems highly unlikely. 

Another explanation might be that the innermost accretion 
disc has different orientation than the larger structures (see, e.g.. 


Navakshinll2005t lUlubav-Siddiki et afll2009t ITremaine & Davisl 

20141) . We can estimate a warped radius to which the accretion 
disc is forced to be in the equatorial plane using a relation by 


ITremaine & Davis! ( 2014 ). which gives 6 x 10 15 cm for the mass 
M ~ 1O 8 M 0 . This is much closer than the distance measured 
for the X-ray absorbing gas (~ xlO 17 cm). This makes the possi¬ 
bility that the innermost accretion disc has different orientation 
than the Broad Line Region. However, a detailed analysis of the 
optical broad spectral lines would be needed for a more relevant 
testing of this hypothesis since warping would have an impact 
on the asymmetry of the broad-line profiles (see, e.g.. lWu et al.l 
120081 ) . 

4.5. Origin of the Broad-Line-Region clouds 

Our detection of an absorbing gas located at the same region 
as BLR may shed light on the origin and geometrical config¬ 
uration of BLR. In general, the broad optical lines are usually 
seen in AGNs viewed at low inclination with a clear view of the 
nucleus. This implies the BLR location close to the equatorial 
plane. However, the origin and exact geometry of BLR is still 
rather unknown and highly debat ed in the literature (see, e.g., 
ISulentic et al.ll2000l:ICraskellll2009L and references therein). 

lElitzur et al.l (12014l) suggested that the BLR and torus are 
most naturall y a continuation of the accretion disc outflow. 
iMiniutti et al.l (12014|) also argued that the BLR clouds might be 
the cold and dense clumps of the smoother warm/hot outflow. In 
Fairall 51, we have also detected highly ionised absorbers with 
an out-flowing velocity comparable to the velocity derived from 
the width of the optical spectral lines. However, the location of 
the highly ionised warm absorber cannot be directly constrained 
from our data since there is no strong evidence for its variability. 
If it is part of BLR, it must constitute a large volume and have a 
very low density. However, it can also originate closer to a black 
hole and be part of the innermost outflow. The higher ionisation 
can be then explained by a stronger illumination from the centre. 

ICzernv & Hrvniewiczl (12011 ) proposed another origin of the 
broad-line clouds. They found that the temperature at BLR is 
close to the dust sublimation temperature, which gives to arise 
dusty winds from the accretion disc. Once they get up from 
the equatorial plane the strong irradiation destroys the dust. The 
BLR clouds lose support from the radiation pressure, which is 
large on dusty particles but much smaller on pure-gas clouds. As 
a result, the clouds are falling back to the disc and again arise. 
This process could explain why the BLR clouds get to the line- 
of-sight of the intermediate type Seyferts like Fairall 51. 

Moreover, the presence of the ionised outflow can provide a 
larger protection of the dusty winds from the destroying irradi¬ 
ation and the clouds may get farther from the equatorial plane. 
Large anti-correlation of the column density of the absorber with 
the intrinsic flux of the nucleus (see Fig.0 is consistent with this 
scenario. Lower flux would allow dusty winds to arise higher 
implying that more clouds may enter into the line-of-sight, and 
thus increase the observed column density of the cold absorber. 
This relation can, however, be also caused by the changes in 
the gas opacity due to the decrease of illumination (see, e.g., 
iKrongold et al.l 120051) . However, with no significant measure¬ 
ments of the change in the ionisation of the absorber, we have no 
quantitative indication for the temperature changes that would 
cause the changes in the gas opacity. 

We note that our location estimation is based on assuming 
the Keplerian orbital velocity. ICzernv et al.l (12015l) showed that 
the vertical velocity of the clouds is of order of z/RxQ k- The z/R 
fraction should correspond to the fraction of sky that is covered 
by the BLR clouds, which is from the obser vations of the order 
of < 0.3 (see, e.g., iGoad & Koratkarlll998 ! ). This satisfies the 
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usage of the Keplerian velocity as the estimated velocity of the 
absorbing cloud. 

Our results also indicate that the distance of the X-ray ab¬ 
sorbing cloud(s) is beyond the warped radius. The warped radius 
may represent an outer radius of the accretion disc, because the 
disc becomes unstable beyond it (see, e.g. JCollin & Zahnl2007l) . 
This allows a possibility that the BLR clouds originate due to a 
fragmentation of still unstable accretion flow. Such clouds would 
not feel strong gravitational attraction towards the disc plane and 
may reach a more inclined orbit. This could be an alternative ex¬ 
planation of the origin of the BLR clouds and why they may get 
into the line-of-sight. 

5. Conclusions 

We performed a monitoring programme of Fairall 51 to con¬ 
strain the characteristic time scale of its X-ray spectral variability 
using the Suzaku satellite. We obtained four new observations 
separated by 1.5, 2 and 5.5 day intervals. A significant spec¬ 
tral change was detected only for the last observation with the 
half flux compared to the first observation, although the flux also 
varied during the individual observations. By comparison of the 
light curves in the X-ray soft (0.5-10keV) and hard (15-55 keV) 
energy bands, we found that the variability occurs almost en¬ 
tirely at the soft X-rays. The 5-days long variability implies the 
location of the absorber to be most likely in the Broad Line Re¬ 
gion. 

Our analysis revealed the spectral complexity in the X-ray 
band. The primary nuclear emission can be well described by a 
power law with the photon index T ~ 2. It is strongly affected 
by at least three absorbers with different ionisations (log £ » 1 
- 4, respectively), the lower ionised one(s) being variable and 
causing the spectral variability. The spectral variability is most 
prominent at 2-6 keV, while the part of the spectrum at < 1 keV 
is probably dominated by a scattered power-law emission com¬ 
ing from either a polar region or an external part of the ionised 
wind itself. 

The hard X-ray spectrum revealed an excess at w 15keV, 
which can be well explained by a reflection from an accretion 
disc that is affected by the relativistic smearing due to strong 
gravitational field around the super-massive black hole and very 
high orbital velocities in the accretion disc. Using the relativistic 
reflection model, we measured the black hole spin, a ~ 0.8, but 
with large uncertainty due to the model complexity. 

The spectral complexity and the rapid absorption variabil¬ 
ity makes Fairall 51a very prominent source for further studies 
with X-ray instruments. In particular, detectors with high-quality 
spectral resolution would help to better constrain the ionisation 
structure of the absorber and its response to the continuum flux 
changes. Future high-resolution observations of Fairall 51 with, 
e.g., planned X-ray mission Astro-H will be therefore desirable. 
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